Finite-time stochastic reduction models 
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A new energy-based stochastic extension of the Schrodinger equation for which the wave function 
collapses after the passage of a fïnite amount of time is proposed. An exact closed-form solution to 
the dynamical equation, vàlid for all finite-dimensional quantum systems, is presented and used to 
verify explicitly that reduction of the state vector to an energy eigenstate occurs. A time-change 
technique is introduced to construct a 'clock' variable that relates the asymptotic and the finite-time 
collapse models by means of a nonlinear transformation. 
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The idea that Brownian motion might play a role in 
models for the collapse of the wave function was fïrst 
envisaged by Wiener and Siegel [1]. Although that spe- 
cific proposal did not in the end fully account for the 
correct probability law, a number of authors have subse- 
quently proposed and developed dynamical collapse mod- 
els driven by Brownian motion that have the key prop- 
erty of being compatible with the probabilistic hypothe- 
ses of Standard quantum mechanics. Progress in this area 
over the last two dècades can be broadly classified into 
two categories: (a) models based on the idea of 'spon- 
taneous localisation' of the wave function [2], and (b) 
models based on the notion of collapse of the state vec- 
tor to an energy eigenstate [3]. See [4] and references 
cited therein for accounts of these approaches. 

The energy-based collapse models, with which we are 
concerned here, have the physical property that the ex- 
pectation of the future energy of the system is given by 
the initial energy. It has been argued on phenomeno- 
logical grounds [5] that the characteristic timescale for 
collapse to an energy eigenstate in such models is of the 
order 
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where AH is the initial energy uncertainty, and it has 
been demonstrated that the choice (1) is consistent with 
empirical observations for a number of different examples 
of quantum systems [6]. The timescale tr is indicative 
of the time it takes for the wave function to reach the 
immediate vicinity of an energy eigenstate. That is to 
say, after the passage of several múltiples of tr, the state 
of the system is, with a high degree of probability, nearly 
indistinguishable from one of the energy eigenstates. It 
should be emphasised, however, that strict collapse, in 
these models, is achieved only asymptotically in time, 
and it has hitherto been unknown how to formulate a 
consistent dynamical collapse model that exhibits a com- 
plete reduction of the state vector in a finite period of 
time. 



The purpose of this article is to introduce a new class of 
energy-based models for which the collapse is completed 
after the passage of a specificd time interval T. Furthcr- 
more, we obtain an exact closed-form expression for the 
state vector process that solves the dynamical equation. 
We show that, remar kably, the finite-time collapse model 
and the Standard infinite-time collapse models are related 
by a time transformation of the form 
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tT 
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where t is the clock time of the finite-time collapse model, 
and t is the clock time in the asymptotic collapse model. 
Thus the finite-time collapse model can be viewed as a 
'fast-forwarded' version of the asymptotic collapse model. 

We begin by stating the main results of the paper. Let 
the dynamics of the state-vector process {\ipt}} be given 
by the following stochastic Schrodinger eguation: 

d\ipt) = -iH\i> t )dt - \a 2 t {È - H t f\^ t )dt 

+ \a t {È - H t )\^ t )àW t . (3) 

Here H t = (iptlHlipt) / (iptl^Pt) is the expectation value of 
the Hamiltonian H in the state \ipt), {W*} is the Wiener 
process, and <j t = aT/(T — t), where a is a parameter. 
Then starting from an arbitrary initial state \ipo) the wave 
function collapses to an energy eigenstate at time T. En- 
ergy is conserved in expectation, and the probability of 
collapse to a state of energy Ei (i — 1,2,3,..., n) is in 
accordance with the Born law 7T; = (^o|IIi|V , o)/(V'o|'0o)> 
where IT is the projection operator onto the Hilbert sub- 
space of states with energy Ei . 

Our model (3) contains two parameters, namely, the 
reduction time T and the energy volatility parameter tr, 
which has the units [energy] -1 [time] -1 / 2 . An application 
of the Ito product rule 



d(X t Y t ) = X t (dY t ) + {dX t )Y t + (dX t )(dF t ) 



(4) 



shows that the normalisation of \ipt) is preserved under 
(3), so d(tjjt\ipt) = 0, and that {H t } satisfies 



1 



where 



V t 



dH t =a t V t dW t , 



(5) 



(6) 



is the energy variance. It follows from (5) that the energy 
process {H t } has no drift and thus satisfies the conser- 
vation law 



E[H U \{W S } 



<8<t\ 



(7) 



for < t < u, where E[— |{W s }o< s <t] denotes conditional 
expectation given the trajectory of the Wiener process 
from time up to time t. 

To vcrify the collapse property we shall examine the 
dynamics of the energy variance process {Vt). In par- 
ticular, by use of the Ito rule (4) together with Jensen's 
inequality we deduce that 



V t < V 



Ja 



VMs, 



(8) 



where V t = E[V t ] is the unconditional expectation of the 
energy variance, which is non-negative. Since Vt <V S for 
s < í, the inequality (8) remains vàlid if we replace V s 
with Vt in the integrand. It follows for any t <T that 



V t <V - er 



tT 



v 2 . 



(9) 



Now suppose Vt > 0. Then there would exist a time 
í = V T/(V +er 2 TV^) < T such that the right side of (9) 
vanishes and hence such that Vt = 0. This contradicts 
our supposition, and hence we conclude that Vr = 0, and 
thus Vt — with probability onc. 

To verify the Born law we must show that the condi- 
tional transition probability 



(10) 



is conserved in expectation, i.e. that ttí — E^] for t < 
T. This follows from an application of Ito's rule, which 
implies that ditn — ert (Ei — H^-iïudWt and thus {tíu} has 
no drift. Because E^t] is the probability of collapse to 
a state with energy Ei, the conservation law shows that 
this is given by the Born probability 7Tj. 

We now procced to investigate the dynamical cquation 
(3) more deeply with a view to obtaining a closed-form 
solution. First we observe that (3) can be cast into an 
integral form that incorporates the initial condition. In 
particular, it follows from (3) that 



\iPt) = Ú t Ml /2 \^) 



(11) 



where Ü t — e lHt is the Standard unitary evolutionary 
operator associated with H, and 



exp (l·l / V a (dW a + a s H s ds) - \H 2 J^ds) 
exp (j*a 8 H a {dW e + a s H s ds) - \ f*er 2 s H*ds) 



(12) 



is a positive operator-valued process. Next we introduce 
a process {W t *} defined by 



W* =W t + I <J s H s ds. 



(13) 



With respect to the Wiener measure P, under which {W t } 
is a Standard Brownian motion, {M 7 ^*} is a Brownian 
motion with drift. Therefore, by the well-known theo- 
rem of Girsanov there exists a measure Q with the prop- 
erty that {Wj*} is a Q-Brownian motion for t € [0, T). 
Since the dynamical law (3) preserves the normalisation 
of \ipt), it follows from (11) that (ipo\Mt\Í>o) — (V'olV'o) 
for t e [0, T), and hence that 



M t = — exp [ H I a^UV- 
where 



\H 2 



a 2 s ds , (14) 



$ t - exp (^J* a s H s dW* s - \ J o 2 s H 2 s d S ^j 

(Vol exp (l·l Jler s dW t *-\H 2 ffódt) |Vo) 



(ipo\ipo) 



• (15) 



We recali that at time t the probability of reduction to 
a state with energy Ei is given by (10). By use of (11) 
and the commutation relation [M t , IT] = it follows then 
that nu — (ipo\ïliMt\ipo) / (tpo\ipo) ■ Substituting (14) into 
this relation we thus obtain 



(VolIT exp (H f*er s dWZ-ÍH 2 J V s 2 d S J |Vo> 
(Volexp (l·l f er s dW s *-±H 2 ffòds) |Vo) 
tt, exp (Ei f* a s dW; - \E 2 f* a 2 ds) 
n exp (Ei jl a s dW* - \E 2 J* a 2 ds) ' 



(16) 



Recalling that er s = erT/(T — s), we find it convenient 
now to introduce a process {Ç t } defined by 



6 = CT 



■dW* 



(17) 



We observe that the right side of (17) is an integral rep- 
resentation for a Brownian bridge [7] in the Q-mcasure. 
Substituting (17) into (16) we infer that the reduction 
probability nu can be expressed in terms of the variable 
£t as follows: 



ití exp ^ ^ — '- — J 



t» exp ( 



<JÍtEíT-±o 2 EïtT 



) 



(18) 



2 



Furthcrmore, on account of the relation H t = ^ i EiiTn 
for the energy process, we deduce that: 



H t = 



^ cx p (- — t±ï — ) 

E* *ï exp ^ T -t J 



(19) 



In view of the exact solution to the asymptotic col- 
lapse model obtained in [8], the fact that H t and nu can 
be expressed as functions of £4 and t suggests there might 
be a simple representation for £ t in terms of elementary 
random data that can be specified without reference to 
the state-vector process. If so, we can substitute such a 
representation for £ t into (19) to obtain a closed-form so- 
lution to the stochastic equation, in place of the integral 
representation (11) which implicitly depends on jJV'í)}- 

We thus proceed as follows: we define a process {fit} 
by the relation fit = £,t — crtHT, where Ht is the ter- 
minal value of the energy. We claim that the random 
variables (3t and Ht are independent. To establish their 
independence it sufficcs to verify that E,[e xl3t+yHT ] = 
E[e xl3t ]E[c yHT ] for arbitrary x,y. Using the tower prop- 
erty of conditional expectation we have 



E[e xfit+yHT ] = E [c xít E \éy- atx ) H T £ 



(20) 



Lct us consider the inner expectation E [ e (y- atx ) H T |£ t j _ 
Using expression (18) for the conditional probability dis- 
tribution of the terminal energy Ht wc obtain 



E 



Jj)—otx)H T 



(21) 



Xy-<rtx)Ei 



exp 



\o 2 E}tT 



Now {$*} is the density process for changing the measure 
from Q to P. That is to say, the expectation E under P, in 
which {W t } is a Standard Brownian motion, is related to 
the expectation E 1 ^ under Q, in which {W t *} is a Standard 
Brownian motion, according to 



E[X u \{W s } <s<t] 



1 



$uX u \{W s } <s<t), (22) 



for any random variable X u that can be expressed as a 
functional of the trajectory {W s } < s <u- Then making 
use of the fact that {£t} is a Brownian bridge under Q, 
we deduce, after some rearrangement of terms, that 



E [. 



= E 7r · 



„ tçr-t) . 



(23) 



Here we have used the facts that if g is a zero-mean Gaus- 
sian random variable with variance 7 2 , then Efe 2 ^] = 
e^ 7 x , and that the variance of the Brownian bridge 
{£t} is t(T — t)/T. This proves the independence of {/?*} 
and Ht- 



The result (23) also establishes that under P the pro- 
cess {Pt} is Gaussian, and has mean zero and variance 
t(T—t)/T for t e [0,T). A similar line of argument shows 
for s < t that the covariance of f3 s and (3 t is 



Cov [f3 s ,[3 t } = s(T-t)/T. 



(24) 



Therefore we conclude that {/3 t } is a P-Brownian bridge. 
As a consequence, there exists a Brownian motion {B t } 
such that {fit} admits an integral representation: 



Pt = (T-t) í -J— dB s . 
Jo T-s 



(25) 



Thus we have verified that the process {£ t } can be ex- 
pressed in the form 



& = atH T + Í3 U 



(26) 



where Ht is the terminal energy, and {(3t} is an indepen- 
dent Brownian bridge. 

Conversely, given a discrete random variable Ht tak- 
ing the vàlues {Ei} with probability {ttí}, and given an 
independent Brownian bridge process {fit} on the interval 
[0, T], we can use the ansatz (26) as a basis for construct- 
ing a closed-form solution to the dynamical equation (3). 

The argument can be sketchcd as follows. First, given 
{£ t } as defined above in terms of H T and {f3 t }, let 7r it = 
¥(Ht = Ei\Ç t ) be the conditional probability that Ht 
takes the value Ei when the value of £ t is specified. It 
follows by use of the Bayes formula that 



F(H T = Ei\Ç t ) = 



irip(Çt\H T = Ei) 
Ei*iP{tt\H T = Ei)' 



(27) 



where 



p(£t\H T = Ei) 



2irt(T-t) y \ 2t{T-t) 



). (28) 



Expression (28) can be deduced from the fact that condi- 
tional on Ht = Ei the variable £t in (26) is normally dis- 
tributed with mean atEi and variance t(T — t)/T. Thus, 
we recover the expression (18) for {ira}- Then by an 
application of Ito calculus we obtain 



d7T ít = a t (Ei - H t )ir lt dW t , 



(29) 



where now the process {Wt} is defined in terms of {Çt} 
by the relation 



W t = J* jj— s (& - o-TH s ) ds + 6 



(30) 



and {H t } is defined by H t = Ei n itEi- The process {W t } 
thus determined turns out to be a Standard Brownian 
motion. This can be demonstrated by showing that {Wt} 
satisfies E[W t |{W u }o<„<J = W s for s < t and {dW t f = 
dt. Finally, if we write 



3 



I0i> = 



n# ) 



(V-oin^o) 1 / 2 



(31) 



for the Standard Lüders states, then the state vector 
{IV't)} that solves (3) can be written in the form 



or, more explicitly, 



-\Ht 



^2vïrïï \<i>i), 



(32) 



where = oT/ (T — s), and £ t = críií T +/3 t . A straight- 
forward calculation shows that \ip t ) — > |</> fc ) as í — > T, if 
we set í?t = E k in (33). By taking the stochastic dif- 
ferential of (33), and using (30), we are led back to the 
starting point (3) , which directly verifies that (33) is the 
solution to the dynamical equation (3). 

We thus conclude that the substitution of (26) in (19) 
yields a closed-form expression for the energy process in 
terms of the exogenously specifiable independent data 
{Ht} and {fit}- As a consequence, we are able to directly 
simulate the solution to the stochastic equation, as well 
as the associated energy evolution {H t }, for an arbitrary 
finite-dimcnsional quantum system. 

Next, we remark on the interpretation of the finitc- 
time collapse model (3). First we observe (cf. [7]) that 
if {B t } is a Standard Brownian motion on the interval 
[0,T], then the process {B T } defined by 



T 



T 



dB s 



(34) 



is a Standard Brownian motion on the time interval 
[0, oo). In particular, let us introduce a new time variable 
by the relation r(t) = tT/(T — t), and define 



T 



'h 



Then the equation (26) can be put into the form 



Vr 



(jtTH T Tf3 t 



T-t T-t 



OTl·ln 



T 



T-s 



dB, 



(35) 



(36) 



on account of the integral representation (25). Making 
use of the relations (34) and t = tT/(t + T), wc thus 
dcducc that 



r] T = cttHt + B T . 



(37) 



However, this is the ansatz (cf. [8]) that determines 
asymptotic collapse to an energy eigenstate. That is to 
say, the expectation H T = E[Ht\t]t] determines the en- 
ergy process {H T } associated with the Standard energy- 
based infinite-time collapse models. Therefore, if we take 



the infinite-time collapse model and replace the time vari- 
able t with a clock r(í), then as t — > T, the clock r(í) 
'ticks' faster and faster. Eventually, in a system for which 
time is measured by 'internal' time t, the collapse takes 
place in a finite time interval T, whereas in a system for 
which time is measured by the clock variable r(í), the 
collapse takes place asymptotically as t — > oo. 

In summary, we have introduced a consistent energy- 
based collapse model that achieves state reduction in fi- 
nite time, thus mceting the 'challenge' posed in Ref. [9]. 
We have verified the collapse property directly by solving 
the dynamical equation (3) explicitly for the state-vector 
process {IV'í}} m terms of independently specifiable ran- 
dom data. We have also obtained closed form expressions 
for the energy process {H t } and the reduction probabil- 
ity process {ira}- An argument based on a time-change 
shows that the finite-time collapse model is related to 
the infinite-time collapse model by a nonlincar transfor- 
mation of the time variable. 
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